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Abstract 

Unstable particles are notorious in perturbative quantum field theory for producing singular propagators in scattering 
amplitudes that require regularization by the hnite width. In this review I discuss the construction of an effective held 
theory for unstable particles, based on the hierarchy of scales between the mass, M, and the width, E, of the unstable 
particle that allows resonant processes to be systematically expanded in powers of the coupling a and E/M, thereby 
providing gauge-invariant approximations at every order. I illustrate the method with the next-to-leading order line- 
shape of a scalar resonance in an abelian gauge-Yukawa model, and results on NLO and dominant NNLO corrections 
to (resonant and non-resonant) pair production of IV-bosons and top quarks. 

Keywords: Unstable particles, effective held theory, perturbative quantum held theory, line-shape, top quark, 
VE-boson 


1. Introduction 

The consistency of the Standard Model (SM) of par¬ 
ticle physics is tested at high-energy colliders primar¬ 
ily through the production and subsequent decay of un¬ 
stable particles. New particles, if discovered, are most 
likely also short-lived. At the current high-energy fron¬ 
tier all known fundamental interactions are perturba- 
tively weak, allowing for very precise theoretical com¬ 
putations in principle. Nevertheless, the application of 
perturbation theory to processes with unstable particles 
is not always straightforward. 

The very notion of an unstable particle requires clar- 
ihcation. In quantum held theory the fundamental en¬ 
tities are the helds from which the Lagrangian is con¬ 
structed, but the excitations of the fundamental helds 
may not correspond to the asymptotic particle states as¬ 
sumed in scattering theory, if they are strongly interact¬ 
ing as is the case for the quarks and gluons of QCD, 
or unstable with respect to decay into lighter particles. 
Relevant cases include the electroweak gauge bosons 
and the top quark, which although all very short-lived, 
have width over mass ratios of a few percent, larger than 
the accuracy of precision calculations. Principal ques¬ 
tions related to held theories with unstable “particles” 


such as their unitarity on the Hilbert space built upon 
the one-particle states of only stable particles have been 
answered many years ago m. The construction of the 
unitary S -matrix is based on certain properties of the ex¬ 
act two-point function of the unstable-particle held. Al¬ 
though a diagrammatic interpretation is assumed, there 
is no explicit reference to a perturbation expansion in 
the coupling that renders the particle unstable. 

Since exact two-point functions are not at hand, this 
raises the question of consistent, successive approxima¬ 
tions. Ordinary perturbation theory in the Lagrangian 
coupling g does not work, since the lowest-order prop¬ 
agator of the unstable particle leads to singularities in 
scattering amplitudes. A well-known remedy of the sin¬ 
gularity is the resummation of self-energy corrections to 
the propagator, which results in the substitution 

p2 _ p2 _ j^2 _ n(p2) ^ ^ 

The self-energy has an imaginary part of order M^g^ ~ 
ME, where E is the on-shell decay width of the reso¬ 
nance, rendering the propagator large but finite. “Dyson 
resummation” sums a subset of singular terms of or¬ 
der {g^M^l[p^ - M^])" ~ 1 (near resonance where 

~ MY) to all orders in the expansion in g^. This pro- 
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cedure leaves open the question of how to identify all 
terms (and only these) required to achieve a specified 
accuracy in and TjM. Failure to address this ques¬ 
tion may lead to a lack of gauge invariance and unitarity 
of the resummed amplitude, since these properties are 
guaranteed only order-by-order in perturbation theory, 
and for the exact amplitude. 

Despite the fact that unstable particle fields have no 
corresponding asymptotic particle states and hence their 
propagators should never be cut, this point is often ig¬ 
nored in practice and the particle is treated in Feynman 
diagram and cross section calculations as if it were sta¬ 
ble (“narrow-width approximation”). This can be justi¬ 
fied when the width is very small, since 


MT 

(p^ - M2)2 -I- M2F2 


7i5{p^ - M^). 


( 2 ) 


The limit holds in the distribution sense and is therefore 
valid only, if the phase-space of the unstable particle is 
integrated sufficiently inclusively, such that the integra¬ 
tion contour in the variable p^ can be deformed far away 
from M^. This is not always the case. An obvious ex¬ 
ample is the line-shape, but also distributions may trap 
the contour near M^. A more accurate treatment than 
Q is also required, when the desired precision exceeds 
the leading-order approximation in F jM. 

Somewhat surprisingly, systematic computational 
schemes to obtain approximations to scattering pro¬ 
cesses involving unstable particles in weak coupling ex¬ 
pansions are relatively recent. The two, which are gen¬ 
eral, are the unstable-particle ejfective theory and the 
complex mass scheme. 


1.1. Unstable-particle ejfective theory 

The singularity of the unstable particle propagator in¬ 
dicates sensitivity to a time scale larger than the Comp¬ 
ton wave-length 1 /M of the particle, evidently its life¬ 
time 1 /F. The presence of two different scales F M 
lies in the very nature of the problem, since a resonance 
with T ~ M would not be identified as such. The main 
idea of unstable-particle effective field theory da is 
to exploit this hierarchy of scales in order to systemati¬ 
cally organize the calculations in a series in the coupling 
g, and F/M. The short-distance scale M is integrated 
out by performing standard perturbative computations 
and the full theory is matched to an effective Lagrangian 
that reproduces the physics at the scale F. The effective 
theory contains a field which describes a resonance 
with momentum p - Mv ■¥ k, where only k ~ F is fluc¬ 
tuating. The resonant field can interact with other soft 
fields with momenta of order F, but off-shell effects at 
the scale M are part of the matching coefficients. The 


power-counting of fields and interactions in the effective 
theory leads to a systematic construction of the expan¬ 
sion in F jM. 

The expansion of amplitudes in matching calcula¬ 
tions is performed around the gauge-invariant location 

iMT (3) 

of the pole in the complex p^ plane corresponding to 
the resonance, where M is identified with the pole mass, 
and F with the on-shell width. The expansion is similar 
to the one performed in the “pole” El El or “double¬ 
pole” (in pair production of resonances) Elll approx¬ 
imation. In a certain sense, unstable-particle effective 
field theory represents the field-theoretic formulation of 
the diagrammatic pole approximation, and generalizes 
it to all orders in perturbation theory and beyond the 
leading power in the F jM expansion. A first step in this 
direction had already been presented in |[8l . 

The effective theory approach is minimal as it iden¬ 
tifies precisely the terms required to achieve a specified 
accuracy in g^, and F/M, and does not include more. 
This makes the calculations particularly simple. Fur¬ 
thermore, the operator interpretation allows for the sum¬ 
mation of large logarithms of F/M through renormaliza¬ 
tion group equations and anomalous dimensions. There 
is a draw-back; the details of the effective theory de¬ 
pend on the inclusiveness of the observable and is not 
valid locally over the entire phase-space, where some 
portions may involve further soft scales of order F in 
addition to {p^ — M^)IM. Even the prediction of the 
resonance line-shape requires matching of the resonant 
(peak) region calculation within the effective theory to 
the off-resonance region computed with standard pertur¬ 
bation theory. 

1.2. The complex-mass scheme 

The complex-mass scheme is an extension of the 
standard on-shell renormalization scheme to unstable 
particles. It defines the complex mass and field renor¬ 
malization constant from the location and residue of the 
pole ([^ of the unstable-particle propagator. The bare 
mass Mo is split into a renormalized mass and countert¬ 
erm through 

Ml ^ Ml + 6Ml, (4) 

and 6MI is part of the interaction Lagrangian and 
treated as a perturbation. The unstable-particle prop¬ 
agator iHp^ - Ml) is never infinite for physical, real 
momenta. The complex-mass scheme was discussed al¬ 
ready in Q, but it was used for the first time in a full 
one-loop calculation only in 2005 IfTOl for the process 
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—» 4 fermions (+y) at high energies, which re¬ 
ceives important contributions from the unstable 
intermediate state. 

Although the standard rules of perturbation theory 
and an expansion in the number of loops apply to the 
complex-mass scheme, a re-ordering and resummation 
of the expansion is implicit, since the propagator is 
of order ll(M^g^) in the resonance region. The assump¬ 
tion is that the complex mass in the propagator captures 
all terms that need to be resummed which is indeed the 
case (see also next section). Since the scheme is only a 
reparameterization of the bare theory, which is not mod¬ 
ified, it is obvious that no double counting occurs. Like¬ 
wise, gauge invariance is assured, since the split Q is 
gauge-invariant and the algebraic identities that guaran¬ 
tee gauge invariance are valid in the presence of com¬ 
plex parameters. Unitarity might be a concern, since the 
unitarity equation involves complex conjugation. How¬ 
ever, since the bare theory is unitary, so must be the 
reparameterized one. What needs to be shown is that 
the theory with the complex-mass prescription is per- 
turbatively unitary in the sense that unitarity violation 
in any given order in the loop expansion are of higher 
order in the expansion parameters (counting YjM ~ g^). 
This point was demonstrated explicitly at one-loop for 
fermion-fermion scattering through a vector-boson res¬ 
onance ifTTIl . and in general in ifT^ . 

The complex-mass scheme is conceptually straight¬ 
forward. It does not require separate treatments of the 
resonance and off-resonance regions, and can easily be 
applied to kinematic distributions. Compared to the ef¬ 
fective field theory method the scheme does not make 
use (explicitly) of expansions in T/M and hence does 
not simplify the problem as much as possible in prin¬ 
ciple. The difficulty of the calculation is equivalent to 
the corresponding standard loop calculation with the 
additional complication of loop integrals with complex 
masses. This is not a practical problem at the one-loop 
order, making the complex-mass scheme the method of 
choice for automated next-to-leading order calculations. 
On the other hand, calculations beyond this order would 
presently be difficult and the resummation of logarithms 
In M/r cannot be performed. 

In the following I do not discuss the complex-mass 
scheme further, but focus on unstable-particle effective 
theory. I use the line-shape of a resonance to illustrate 
the framework and the discuss results on pair produc¬ 
tion of W-bosons and top quarks near threshold which 
(I believe) benefit particularly from this method. 


2. Line-shape of an unstable particle 

In this section, which follows we consider a 

toy model that involves a massive scalar field, cp, and 
two fermion fields. The scalar as well as one of the 
fermion fields, ijj, (the “electron”) are charged under an 
abelian gauge symmetry, whereas the other fermion, x, 
(the “neutrino”) is neutral. The model allows for the 
scalar to decay into an electron-neutrino pair through a 
Yukawa interaction. The model describes the essential 
features of the Z-boson line-shape in the SM llT3]l . Its 
Lagrangian is 

-£ = (Df^tpy D^tp - M^(f>^tpXifix 

+yHx + y'’4>^x^ - + -Cct, (5) 

where M denotes the renormalized mass, not neces¬ 
sarily the pole mass M defined by Q, Xct the coun¬ 
terterm Lagrangian, and - igA^. We define 

ag = g^l{An), a,, = (yy*)l{An) (at the scale yu) and as¬ 
sume ag ~ Oy ~ a, and = 4/(47r) ~ a^/(47r). 

The line-shape is the totally inclusive cross section 
for the process 

v{q) + e^{p) X (6) 

as a function of s = (p -i- q)^, which can be calculated 
from the imaginary part of the forward scattering ampli¬ 
tude 7”(s){^In particular, we are interested in the region 
s ss M^, or more precisely s - ~ MY ~ aM^ « M^, 

where we expect an enhancement of the cross section 
due to the resonant production of the scalar. Defining 
the dimensionless variable 


the cross section far away from the resonance can be 
expanded in g^ in the usual manner according to 

cr^g^M6) + g^f2(6) + .... ( 8 ) 

At every order, the coefficient f„(6) is a function of the 
variable 6 ~ 1. On the other hand, near resonance we 
may exploit 5 1 to expand the amplitude in 6. At the 

same time, as g^ jd ~ 1 since T ~ Mg^, some terms must 
be summed to all orders. A systematic approximation 


*The total cross section of process l|^ is not infrared finite for 
massless electrons due to an initial-state collinear singularity, which 
has to be absorbed into the electron distribution function. In what 
follows it is understood that this singularity is subtracted minimally. 
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to the line-shape in the resonance region therefore takes 
the form 

( 2 x" 

^ x{l(LO);g2,5(NLO),...) 

' 

= + + (9) 

with non-trivial functions hn(g^ld) at every order in the 
reorganized expansion. The effective theory identifies 
the relevant terms and constructs the expansion 0. 




Figure 1: Reduced diagram topologies in 2 —» 2 scattering near 
resonance. Left: resonant scattering. Right: non-resonant scattering. 


2.7. Relevant modes and reduced scattering diagrams 

The effective theory is based on the hierarchy of 
scales r M. In a first step we integrate out hard mo¬ 
menta k ~ M. The effective theory will then not contain 
any longer dynamical hard modes since their effect is 
included in the coefficients of the operators. The hard 
effects are associated with what is usually called fac- 
torizable corrections, whereas the effects of the dynam¬ 
ical modes correspond to the non-factorizable correc¬ 
tions ISl. On the level of Feynman diagrams, the hard 
contribution can be identihed directly using the method 
of regions to separate loop integrals into various contri¬ 
butions na. The hard part is obtained by expanding the 
full-theory integrand in 6. 

The modes to be described by the effective La- 
grangian correspond to kinematically allowed scatter¬ 
ing processes with virtualities much smaller than M^. 
Particles with masses above MY are no longer present, 
except for the unstable particle, which by construction 
is close to mass-shell. To account for this, we write 
the momentum of the scalar particle as P = Mv -t- k, 
where the velocity vector v satisfies = 1 and the 
residual momentum k scales as M6 ~ T. In analogy 
to heavy-quark effective theory (HQET) we remove the 
rapid spatial variation from the (f> field and define 

(pv{x) = e'^'’’^ V+tpix), where V+ projects onto the pos¬ 
itive frequency part to ensure that is a pure destruc¬ 
tion field. A held with momentum huctuations k ~ Y is 
called a “soft” held. Thus, for the soft scalar held (p^, we 
have ~ M^5. This remains true if the scalar 

particle interacts with a soft gauge boson with momen¬ 
tum M6, so the effective Lagrangian should contain soft 
(s) helds for every massless held of the full theory. 

The unstable particle is produced in the scattering of 
on-shell particles with large energy of order M. These 
can remain near mass-shell by radiating further ener¬ 
getic particles in their direction of Hight. The effective 
Lagrangian must therefore also contain hard-collinear 
(cl) modes with momentum scaling 

n+p ~ M, Pj_~ n_p ~ MS (10) 


for all massless helds of the original Lagrangian. Here 
are two light-like vectors with n+ ■ n_ = 2, n_ is 
the direction of the electron four-momentum, and is 
transverse to n_ and n+0 

The space-time picture of the kinematically allowed 
processes is very simple and the corresponding re¬ 
duced diagram topologies are shown in Figure [T] for 
the forward-scattering amplitude. The left diagram de¬ 
scribes the production of the resonance through a hard 
process, represented in the effective theory by some lo¬ 
cal operator (9®, and its subsequent propagation over 
distances of order 1 /F. The resonance (double line) can 
interact with soft huctuations. The initial-state electron 
leg can be dressed with collinear corrections. However, 
collinear modes cannot be exchanged across the dou¬ 
ble line, since this would not leave enough energy to 
produce the scalar near resonance. The process just de¬ 
scribed is represented in the effective theory by the hrst 
line of the master formula 


iT = 


kl ^ 


ctx <ve| T(iOf(0)iOf(x))\ve) 


<ve|/<9®(0)|ve). 


( 11 ) 


for the forward-scattering amplitude. 

The scattering may also occur without the produc¬ 
tion of the scalar near its mass-shell (right diagram in 
Figure [TJ. In the present toy theory this still requires 
an intermediate scalar line, since the neutrino has only 
Yukawa interactions. The scalar may be off-shell, be¬ 
cause the electron has radiated an energetic (hard or 


^In the general case several types of collinear modes are required, 
one for each direction defined by energetic particles in the initial and 
final state. For the inclusive line-shape we calculate the forward¬ 
scattering amplitude, so no direction is distinguished in the final state. 
We then need two sets of collinear modes, one for the direction of the 
incoming electron, labelled by “cl” (or often simply “c”), the other 
for the direction of the incoming neutrino (labelled “c2”). Since the 
neutrino is electrically neutral, the collinear fields if/c 2 , Ac 2 and^ci 
appear only in highly suppressed terms, so we can ignore them here. 
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collinear) photon before it hits the neutrino. In this case 
the invariant mass of the colliding electron-neutrino sys¬ 
tem is of order but not near producing a non¬ 
resonant scalar. In the effective theory this process is 
represented by a local four-fermion operator with¬ 
out 01, fields. In general, non-resonant scattering in¬ 
cludes all “background processes”, which produce one 
of the final states under consideration. This topology 
does not involve a resonant heavy scalar, and both soft 
and collinear fields can be exchanged across the dia¬ 
gram. The matrix elements in ( [TT] i are understood to 
be evaluated with the effective Lagrangian. 

2.2. Construction of the ejfective Lagrangian 

We divide the effective Lagrangian into three parts. 
Roughly speaking, the first, -Chset. describes the heavy 
scalar field near mass-shell and its interaction with the 
gauge field. The second part, X±, describes energetic 
fermions and their interactions with the gauge field. Fi¬ 
nally, the third part, -Cint contains the local operators 
and responsible for the production of the resonance 
and off-shell processes. In the following, we write down 
all terms needed for a next-to-leading order (NLO) cal¬ 
culation of the line-shape. 

The soft Lagrangian -Chset is an extension of the 
HQET Lagrangian mi to a (here scalar) particle whose 
mass-shell is defined by the complex pole location Q. 
The residual mass term which is usually set to zero in 
HQET by choosing M to be the pole mass of the heavy 
quark, is now necessarily non-vanishing and complex. 
The relevant terms are 


-Chset - 2M(f>l (iv ■ Ds - — 


■2Mfl 


(iDs,Tf 

2M 


Ad) 

[Ad)]2 

8M 


A(2) 


- + fsiPsfs +Xsi ?Xs, 


( 12 ) 


where ifis (Xs) denotes the soft electron (neutrino) field 
and the covariant derivative 0^ = 8- igA^ includes the 
soft photon field. Eurthermore, a^ = a^ - {v ■ a) for 
any vector. The only non-trivial short-distance match¬ 
ing coefficients in this expression are the quantities 
to be defined below. 

The bilinear terms in the soft scalar field fy are de¬ 
termined by the requirement that Xhset reproduces the 
two-point function of the scalar in the full theory close 
to resonance. Denoting the complex pole of the propa¬ 
gator by and the residue at the pole by R^, the prop¬ 


agator near resonance can be written as 


iR^ _ iR,!, 

P^-Ml ~ 2Mv -k + k^-{Ml- m 2) ■ 


(13) 


We now define the matching coefficient A s {Ml — 
M^)IM. There are two solutions to - Ml, one of 
which is irrelevant since it scales as v ■ k ~ M. Eor the 
other we find 


vk 


-M + 



-I- MA - k\ 


A 

2 


A^ -H 44 
8M 


+ C)(d3), 


(14) 


where we expanded in 6 in the second line, using A ~ 
k-y ~ M6. Expanding A = A*-'^ into terms of order 

g^‘, we deduce the bilinear terms in (12i from the dis¬ 


persion relation (141. Gauge invariance of the effective 
Lagrangian implies that the leading soft-photon inter¬ 
actions can be obtained from the bilinear terms by re¬ 
placing 8 —» Dg. The gauge invariance of the matching 
coefficient follows from the invariance of the unstable- 
particle pole M*. 

In the underlying theory the full renormalized prop¬ 
agator of the unstable particle is given by i{s - M^ - 
n(i)) ', where -/n(i) corresponds to the amputated 
IPI graphs including counterterms. Comparing this to 
(13 I and expanding n(i) around M^ and in the number 
of loops in the form n(s) = M^ Yjk.i H*^*’^*, where it is 
understood that ~ we obtain 

A = MH'*’®) H- M(n<2’°) -H -H .... (15) 


n*''") and n<2.0) + (but not n(2.0) and H*'-') 

separately) are infrared-finite, which justifies the inter¬ 
pretation of A as a short-distance coefficient. Explicit 
results for A*^'' and A*^^' in the MS and pole renormal¬ 
ization scheme can be found in ||3l- Here we only note 
that in the pole scheme (M s M), we have A = -/E, in 
which case the residual “mass” is purely imaginary and 
coincides with the on-shell width. 

Each term in -Chset c™ be assigned a scaling power 
in 6. Since D, ~ k ~ T ~ M6 and A^b ^^2 ^ 
M6, both terms in the first line of ( [T 2 I 1 are of equal size 
and leading terms. The unstable-particle propagator is 
therefore 


2M(vk-A(b/2)’ 

which corresponds to a fixed-width prescription. The 
linearity of the propagator in the (residual) momentum 
makes calculations in the effective theory particularly 
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simple. The fact that only appears in the leading- 
order Lagrangian proves that only the two-point func¬ 
tion in the original theory needs to be resummed by in¬ 
cluding the one-loop self-energy into the unperturbed 
Lagrangian. No higher-point functions require resum¬ 
mation, which is intuitively obvious, since the origin 
of the long-distance scale is associated with a single¬ 
particle effect, the life-time of the resonance. 

In momentum space the propagator ( [Thl l of the (fiy 
field scales as l/d. Hence, because J d'^k counts as 6^, 
the soft scalar field (p^ix) scales as It follows that 
the terms in the second line of ( [l^ scale as 6^. Being 
suppressed by one power in 6 or^ relative to the first 
line, they must be included only in a calculation of the 
line-shape with NLO precision. Finally, since scales 
as 6 and the soft fermion fields scale as 6^^^, the terms 
in the last line of (12 1 scale as and represent leading 
interactions among the soft, massless modes. By adding 
further terms the Lagrangian can be improved to any ac¬ 
curacy desired. 

Next, we turn to the construction of the effective La¬ 
grangian, X.±, associated with the energetic fermions. 
The interactions of collinear modes with themselves and 
with soft modes are described within soft-collinear ef¬ 
fective theory (SCET) |[T6l [HI HS] [H . The coupling 
of collinear modes to the scalar field cp^, and among 
collinear fields with different directions produces off- 
shell fluctuations, which are not part of the effective 
Lagrangian. The momenta associated with generic 
collinear fields tp^i andT'c 2 do not add up to a momen¬ 
tum of the form P - Mv -H k. This kinematic con¬ 
straint is implemented by adding the production and 
non-resonant operators, (9® and respectively, as 
external “sources” for the specific process. The line- 
shape is then given by the correlation function ( 111 . 

Alternatively, the dynamical hard-collinear modes 
can be integrated out in a second matching step, in 
which the collinear functions (labelled “C” in Figure [^l 
appear as matching coefficients of (non-local) operators. 
The new effective Lagrangian contains an “external- 
collinear” electron mode with momentum Mn_/2 -i- k, 
which describes the remaining soft fluctuations k ~ 6 
around the fixed large component. Similar to the reso¬ 
nance field, we extract the large component and define 
ifinJx) = e‘i4l2(n-x)rp^ ^c\{x), where projects on the 
positive frequency part of \pc\- Adding the correspond¬ 
ing field with n_ and n+ exchanged for the neutrino, the 
soft interactions of the external-collinear field are given 
by 


-C± = ifiniri^D,— [fin 


+ — Xn^- 


(17) 




Figure 2: Scalar self-energy correction to the forward-scattering am¬ 
plitude. 


at leading power. 

With the external-collinear modes we can implement 
the production and non-resonant sources as interaction 
terms in Xjnt. At NLO the relevant terms read 

-Cint = Cy(pylpn_Xn^ + Cy*^lxn^lj/n- 

+ {ipn^XnJiXny ^PnJ, ( 18 ) 

where C = 1 and D are the matching coefficients. 

The two lines correspond to the two reduced diagram 
topologies in Figure We note that the effective La¬ 
grangian is not manifestly hermitian, since it describes 
the decay of the scalar. Nevertheless, it generates a uni¬ 
tary time evolution, since it reproduces by construction 
the unitary underlying theory to the specified order in 
the expansion in 6. 

The external fields scale as Thus, an insertion 
of a 4>4iX operator results in J d^x<pyip„ Xn+ ~ The 
forward-scattering amplitude requires two insertions of 
this operator. Accounting for the scaling of the exter¬ 
nal state {ve~\ ~ 5'', we find ~ jd for the am¬ 
plitude at leading order, which is the expected result. 
The four-fermion operator is suppressed in 5 and results 
in a contribution of order g^ to T. Thus, to compute 
the NLO correction 7“*'* we need C*'', the 0{g^) con¬ 
tribution to the matching coefficient C, while D is only 
needed at tree level. The matching coefficients are ob¬ 
tained from the hard contributions to the corresponding 
on-shell three- and four-point functions in the full the¬ 
ory. I refer to El for the precise matching equation as 
well as the explicit results. 

2.3. Example diagram 

It is instructive to discuss how the self-energy cor¬ 
rection to the intermediate scalar in the full theory, see 
Figure]^ is represented in the effective description. We 
first separate the hard and soft contributions the one- 
loop self-energy, n(i) = n/,(i) -i- 0^(5), and then expand 
the hard part n/,(i) = The soft part is 

reproduced by the effective theory self-energy. The first 
term in the hard expansion is gauge-invariant and 
contributes to A*^*^ as discussed before. This term is al¬ 
ready relevant to the leading-order line-shape. The next 
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term cancels one of the adjacent scalar propaga¬ 
tors, such that the self-energy correction merges with 
the local production vertex, is gauge-dependent. 

The gauge dependence cancels with the vertex diagram 
to produce a gauge-independent NLO hard-matching 
coefficient Continuing in this way, we find that 

contributes to the one-loop matching coefficient 
of the four-fermion operator {^n-Xn+){Xn^^n-), be¬ 
cause the scalar propagators to the left and right are both 
cancelled. The contribution is again required to ob¬ 
tained a gauge-invariant one-loop matching coefficient 
a, though it is already a NNLO term for the line-shape. 

This example illustrates the power of the effective 
field theory method. It automatically breaks a diagram 
into different pieces and organizes them into gauge- 
invariant objects. The power-counting associated with 
the Lagrangian allows one to identify the terms relevant 
for a specified accuracy before any explicit calculation 
needs to be performed. 

2.4. Line-shape at next-to-leading order 

Our goal is to carry out this programme for the 
forward-scattering amplitude at NLO, where 

7“*®^ sums up all terms that scale as ~ 1 and 

7“*'* contains all terms that are suppressed by an addi¬ 
tional power of g^ or 6. At leading order there is only 
one diagram, involving two three-point vertices and one 
resonant scalar propagator. We get 

* 7 "® = [Hp)v(q)] MqMp)], ( 19 ) 

2MD 

where we defined D = sfs - M - A*^'V2. The inclusive 
line-shape is related to 7“^®^ by cr = ImT^^^/s through 
the optical theorem. The above expression gives a Breit- 
Wigner distribution in yfs. 

In the effective theory there are three classes of 
diagrams that contribute to corresponding to 

hard, hard-collinear and soft contributions. The hard- 
collinear corrections to the external lines lead to scale¬ 
less integrals and vanish. The hard corrections consist 
of a propagator insertion aproduction 

vertex insertion and a four-point vertex diagram 
due to the operator in Xint, as shown in the up¬ 

per diagrams of Figure The sum of these diagrams 
reads 




2C<b- 


[A(1)]2 ^ ^(2) ^ 

8DM 2D 2M 


( 20 ) 


The soft-photon one-loop corrections (lower set of di¬ 
agrams in Figure computed in the effective theory 




Figure 3: Hard (upper) and soft (lower diagrams) contributions to 


result in 


iT 


( 1 ) _ ;rri0) 


X 


(4;r)2 


- 4L H- 


5n^ 


( 21 ) 


with L — ln(-2£)//i). The partonic line-shape is ob¬ 
tained after subtracting the initial-state collinear singu¬ 
larity and taking the imaginary part. The partonic line- 
shape must then be convoluted with the electron distri¬ 
bution function. 

We note the simplicity of the result, which is a conse¬ 
quence of the fact that the complete calculation is bro¬ 
ken into separate single-scale calculations by factoriz¬ 
ing the hard and soft regions. The NLO correction leads 
to a distortion of the line-shape relative to the Breit- 
Wigner form, which in non-inclusive situations can de¬ 
pend on the final state. Fitting a measured line-shape 
to the Breit-Wigner form rather than the true shape pre¬ 
dicted by theoretical calculations leads to errors in mass 
determinations. In the present toy model, choosing the 
pole mass M - 100 GeV (such that the MS mass is 
M = 98.8 GeV at LO and M = 99.1 GeV at NLO) and 
couplings g^l{4n) = \y'^ jiAn) - 0.1 to mimick the pa¬ 
rameters of electroweak gauge bosons, the error would 
be of order 100 MeV. 

Figure]^ shows the leading-order partonic line-shape 
in the effective theory and the tree-level (order a^) cross 
section off resonance in the full theory. The two re¬ 
sults agree in an intermediate region where both calcu¬ 
lations are valid. This allows to obtain a consistent LO 
result for all values of xfs. The figure also shows the 
NLO line-shape for the numerical values given above. 
In order to obtain an improved NLO result in the en¬ 
tire region of sfs, the NLO line-shape would have to be 
matched to the NLO off-resonance cross section in the 
full theory. 

The method discussed here makes NNLO line-shape 
calculations in 2 ^ 2 scattering possible with present 
techniques. An outline of such a calculation has been 
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Figure 4: The line-shape (in GeV“^) in the effective theory at LO 
(light grey/magenta dashed) and NLO (light grey/magenta) and the 
LO cross section off resonance in the full theory (dark grey/blue 
dashed) as a function of the centre-of-mass energy (in GeV). Figure 
from (3- 


given in |[3l, though no complete calculation has been 
performed to date. 

3. Pair production near threshold 

I reviewed in some detail the case of the line-shape, 
since it serves well to illustrate the general framework 
of unstable-particle effective theory. However, some of 
the more interesting results concern pair production of 
unstable particles, specifically the W bosons and top 
quarks, near threshold. In collisions very precise 
measurement of the masses of these particles can be ob¬ 
tained from a threshold scan. 

The threshold dynamics is determined by the inter¬ 
play of the strength of the electromagnetic (W bosons) 
or colour (top quarks) Coulomb force and the size of the 
decay width of the particle. The small parameters are 

6 =^, v^ = {^fs-[2M + /r])/M, (22) 

and the coupling a - HAn). For W bosons, Fiv ~ 
Mn/ttEw and therefore the effective strength of the 
Coulomb force is a^mlv ~ ^ 1. This leads to an en¬ 

hancement, but the Coulomb force is never (9(1), and no 
resummation is needed Il20ll . The rapid decay of the W 
boson prevents the formation of any visible reso¬ 

nance. The situation is different for top quarks, since the 


Coulomb force is generated by QCD, while the decay 
still occurs through the electroweak interaction. Count¬ 
ing as ~ we now find a^/v ~ 1. Diagrammati- 
cally, ladder diagrams that contain these enhanced terms 
must be summed to all orders in perturbation theory, 
which generates toponium bound-states in the spectral 
functions. Since the characteristic energy near thresh¬ 
old E ~ Mv^ is of order T, the bound-states appear as 
broad resonances, of which only the first one leaves a 
distinctive feature in the tt cross section MM- 

In the following I review results for W and top pair 
production near threshold obtained within the effective 
field theory approach, leaving out all the technical de¬ 
tails that can be found in the original papers. 


3.1. W bosons 

This subsection summarizes results from Il23l l24l . 
We consider the process e~e^ —> jjTVnudX with centre- 
of-mass energy sfs - 160... 170 GeV, where it is domi¬ 
nated by a W^W~ intermediate state near threshold with 
subsequent semi-hadronic decay. The inclusive cross 
section is extracted from specific cuts of the forward 
amplitude 

d- = ^ ImV[(e‘e'" ^ (23) 

which also includes diagrams with only a single internal 
W line. We perform a “QCD-style” calculation of the 
“partonic” cross section cr with massless electrons in 
the MS scheme, and convolute it with the MS electron 
distribution function; 


cr(s) = I dXidX2 fe/e(x\) fe/eixi) d-(XiX2S). (24) 

Jo 


The MS electron distribution function depends on m^, 
but not on sfs, M, F. 

In the effective field theory (EFT) the W bosons are 
described by two non-relativistic three-vector fields QJj, 
where a - ± refers to the charge of the W. The HSET 
Lagrangian relevant to a single (scalar) unstable parti¬ 
cle is replaced by the PNRQED Lagrangian 1251, gen¬ 
eralized to the case of an unstable vector particle. The 
relevant terms are 


-Cpnrqed - 2^ 




iD° + 




2M, 




Q' i 

ml “ 


-J(Pr [nl'QL](.x -H r) [G+G{](.r). (25) 
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cr[fbj 



Figure 5: Successive EFT approximations: LO (long-dashed/blue), 
(dash-dotted/red) and NLO (short-dashed/green). The 
solid/black curve is the full Born result computed with Whizai'd/ 
CompHep. The EFT approximation is indistinguishable from 

the full Bom result on the scale of this plot. Figure from ED 


The master formula for the forward amplitude ^ co¬ 
incides with but the production and non-resonant 
operators are now of the form 

Of = Cf , (26) 

of = cf (e,,rie,y(e,p2ec,), (27) 


with Ti, r 2 Dirac matrices, = a'bj -i- a^b', and n 
the unit-vector in the direction of the incoming electron 
three-momentum. 

Due to the 1/v enhancement of electromagnetic 
Coulomb exchange, the systematic expansion of JTI goes 
in powers of Vd. Also, the non-resonant term ap¬ 
pears as such a “N^^^LO” correction, since the lead¬ 
ing imaging parts of Cf are proportional to a^, while 

Vdin 

The EFT constructs an expansion in T jM and { yfs — 
2M)IM of the full theory Born cross section. Before 
turning to radiative corrections it is instructive to com¬ 
pare successive terms in this expansion to the full Born 
result computed numerically (using Whizard ll26l and 
CompHep lIZTl '). This is shown in Figure The FO 
non-relativistic approximation overestimates the true re¬ 
sult. The N'/^lO non-resonant correction yields a 
(nearly) constant, negative term and provides already 
good agreement close to the nominal threshold at a/s « 
161 GeV. To extend the approximation in a wider re¬ 
gion around the threshold, it is necessary to include all 
terms up to N^^^FO. 



Figure 6: One-loop diagram for the hard-matching coefficient (left) 
and a soft NLO contribution to the forward-scattering amplitude in 
the elfective theory (right). 


While constructing an expansion of the Born cross 
section when an exact, numerical result is readily avail¬ 
able, appears as an unnecessary complication, the com¬ 
putation of the NFO radiative correction in unstable- 
particle effective theory 12^ is remarkable simple com¬ 
pared to the corresponding calculation in the complex- 
mass scheme IIT0ll28ll . The most complicated part is the 
computation of the NFO matching coefficient of the op¬ 
erator (261, which, however, is a standard one-loop cal¬ 
culation. A representative diagram is shown in Figure]^ 
left. The diagram on the right displays a soft, “non- 
factorizable” NFO correction to the two-point function 
of production operators in (fTTll, and results again in a 
simple expression, similar topTb. A comparison with 
the complex-mass scheme calculation and the double¬ 
pole approximation (DPA), including QCD corrections 
and initial-state radiation is given in the following ta- 
ble|3 The numerical difference of 1% between the EFT 
and full ee4f results is presumably in part due to the 
N^^^FO correction associated with the NFO matching 
coefficient of Of, which is implicitly contained in the 
NFO full ee4f calculation. 



cr(e e'^ p v^MJX)(fb) 

Vi [GeV] 

Born (SM) 

EFT 

full ee4f 

DPA 

161 

107.06(4) 

117.38(4) 

118.77(8) 

115.48(7) 

170 

381.0(2) 

399.9(2) 

404.5(2) 

401.8(2) 


We can now estimate the theoretical uncertainty in 
the W mass determination from a threshold scan. Fig¬ 
ure shows K = cr{s,Mw + 6Mw)lcr{s, Mw) for 
Mw - 80.377 GeV and different values of 6Mw as 
function of the cms energy. The relative change in 
the cross section is shown as dashed lines for 6Mw - 
+ 15,+30,+45 MeV. The shape of these curves shows 
that the sensitivity of the cross section to the W mass is 
largest around the nominal threshold Vi 161 GeV, as 
expected, and rapidly decreases for larger (The loss 
in sensitivity is partially compensated by a larger cross 
section, implying smaller statistical errors of the antici¬ 
pated experimental data.) The shaded areas provide an 


^The factor Vt arises from the leading-order EFT matrix element 
and corresponds to the phase-space suppression near threshold. 


'^The “full ee4f” column refers to the erratum of Ei 
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K 



VF[GeV] 


Figure 7: W-mass dependence of the total cross section. All the cross 
sections are normalized to (t{s, M]/\/ = 80.377 GeV). See text for ex¬ 
planations. Figure from E). 


estimate of the uncertainty from uncalculated 
terms. The inner band is associated with the interfer¬ 
ence of single Coulomb exchange with one-loop hard 
or soft corrections, which are genuine NNLO correc¬ 
tions in other schemes. The outer band accounts for the 
non-resonant term already mentioned above. Finally, 
the line marked “ISR” estimates ambiguities in the im¬ 
plementation of initial-state radiation. This represents 
the largest current uncertainty. In order to obtain a com¬ 
petitive determination of M^, one eventually needs a 
more accurate computation of the electron distribution 
function. 

Since this is not a fundamental problem and since 
the full theory NLO ee4f calculation is available, the 
accuracy of the theoretical prediction is limited by the 
terms in the 6 expansion, which correspond 
to two-loop corrections (in the complex-mass scheme). 
Some of the diagrams together with their EFT represen¬ 
tation are shown in Figure These consist of mixed 
hard-Coulomb corrections (first column), interference 
of Coulomb exchange with soft and collinear radiative 
corrections (2nd and 3rd column, respectively), and a 
correction to the electromagnetic Coulomb potential it¬ 
self. These genuine higher-order corrections have been 
computed ll^ and were found to be below 0.5%, lead¬ 
ing to shifts of W mass of less than 5 MeV. 

Up to now, we considered the total cross section for 
the flavour-specific final state p Since experi¬ 

mentally certain cuts must be applied, it would be desir¬ 
able to compute directly the cut cross section in the EFT. 
A framework to implement arbitrary cuts while main- 



Figure 8: Some NNLO diagrams that count as N^^^LO in the 6 expan¬ 
sion and their EFT representation (upper line). 


o-lfb] 



Figure 9: Comparison of the Born cross section in the full SM at = 
161 GeV computed with WFIIZARD (red dots) with the effective- 
theory result for the loose-cut implementation (dashed blue curve) and 
the tight-cut implementation (solid black curve). Figure from (Ml 


raining an expansion in the power-counting parameter 6 
is not available and probably difficult to achieve. The 
specific case of invariant-mass cuts - M^| < 
on the VT-decay products has been considered in ll24l . 
The implementation depends on how A scales with the 
parameter 6. For loose cuts, A ~ Since by as¬ 

sumption the virtualities in the EFT are at most of or¬ 
der MY ~ M6, the loose cut does not affect the EFT 
diagrams. However, the hard-matching coefficients are 
modified and acquire a dependence on A in addition to 
the other short-distance scales. The situation is reversed 
for tight cuts with A ~ MY ~ M V^. The tight cut cuts 
into the (approximate) Breit-Wigner distribution of the 
single-W invariant mass distribution and therefore must 
be applied to the calculation of the EFT loop integrals. 
On the other hand, it eliminates off-shell contributions, 
and hence the short-distance coefficient C® of the non¬ 
resonant four-electron operator (27 1 vanishes. Figure]^ 
shows good agreement of the effective-theory calcula¬ 
tion of the cut Born cross section with the numerical 
result from WHIZARD in the regions where the respec¬ 
tive loose/tight-cut counting rule is appropriate. 
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Figure 10: Top quark pair production cross section near threshold, 
normalized to Ana\^l{?)s) at NNLO in non-relativistic, resummed 
perturbation theory in the PS mass scheme i29] . The width of the 
band reflects the theoretical uncertainty estimated from scale variation 
in the indicated range. m( Ps(20GeV) = 171.5 GeV, T, = 1.33 GeV. 


3.2. Top quarks 

The strong Coulomb attraction of the coloured top 
quarks, when their relative velocity is small near thresh¬ 
old, requires the resummation of certain QCD correc¬ 
tions to all orders in perturbation theory. A system¬ 
atic formalism employs a sequence of matching steps 
to define hard matching coefficients and potentials in 
non-relativistic effective field theory. The ingredients 
required up to the third order in non-relativistic pertur¬ 
bation theory are reviewed in ISO). Figureshows the 
second-order (NNLO) result from BTIl . which exhibits a 
toponium resonance slightly below the nominal thresh¬ 
old. Note the order counting here is such that v ~ 


defines the expansion parameter. The top-quark width 
r,/m, ~ ffiEW ~ a] second-order in this counting. 

QCD predictions of the top-pair production cross sec¬ 
tion such as the one shown in the figure are based on the 
calculation of QCD correlation functions with the sub¬ 
stitution E - ^[s - 2m, E + iF, to account for the 
top-quark width 112111221 . This prescription corresponds 
to computing the first (resonant) term in (11 1 with 
given by the non-relativistic top-quark (axial-) vector 
current, and with an effective Lagrangian that accounts 
for the width through A = -iF, in the first bilinear term 
in ( [25| ) (adapted to quarks), but not in the further kinetic 
corrections. 

The limitations of this approximation manifest them¬ 
selves within the (NR)QCD calculation itself. The cur¬ 
rent correlation function G{E) exhibits an uncancelled 
ultraviolet divergence from an overall divergence of the 
form [dG(£’)]overaii (XsEje in dimensional regulariza¬ 
tion {d - A - 2e) 1321. Since E acquires an imaginary 
part r, ~ m,aEw, the divergence survives in the cross 
section. 


Im [^G(£)]overall IJlf X (28) 

e 

and appears first at NNLO (since at LO, G{E) ~ v ~ a^). 
A consistent calculation therefore requires that one con¬ 
siders the process —> W^W^bh within unstable- 

particle effective theory including the effects of off-shell 
top quarks and processes that produce the W^W^bh fi¬ 
nal state with no or only one intermediate top-quark 
line. The two terms of ([TT]i can be identified with 


^e+e-^W+W-bb - O'e*e - + (^e*e--^W*W-bbn„„J.rw}- (29) 

pure (NR)QCD 


Both terms separately have a “finite-width scale depen¬ 
dence” related to the uncancelled 1/e poles, and only 
the sum is well-defined. For consistency, both terms 
have to be defined with the same (dimensional) regular¬ 
ization prescription. 

While the explicit finite-width scale dependence is 
seen first at NNLO, the leading non-resonant contribu¬ 
tion already appears at NLO. Somewhat surprisingly, 
this was realized only recently. At this order the match¬ 
ing coefficient of the non-resonant operator is equiva¬ 
lent to the dimensionally regulated pro¬ 

cess with F, = 0, expanded in the hard region around 


s - Amj. The corresponding calculation has been 
performed in two independent ways [33] El- In ll^ 
invariant-mass cuts m, - AM, < M,j < m, + AM, on the 
decay products of the (anti-) top quark can be included 
following the method discussed above for W bosons. 

It is convenient to represent the calculation of the cut 
two-loop diagrams contributing to —> bW^t in the 

form 

r ' dp^ im) - p2) ¥ ^^ j (30) 

Ja2 'm, m, / 

leaving the integration over p^ = {pt, + pw*Y to the end. 
The lower limit depends on the invariant-mass cut and 
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Figure 11: Hai'd, off-shell top diagram contribution to e^e —> bW^t, 
which leads to a linear infrared divergence. 



Figure 12: Relative sizes of the non-resonant corrections with re¬ 
spect to the tt LO cross section in percent: NNLO singular terms 
(upper blue lines) and NLO crl^on-TcJ (lower black 
lines). Solid (dashed) lines correspond to an invariant-mass cut AMr = 
35 GeV (AM[ = 15 GeV). Figure from 1351 . Here the top pole mass 
m{ = 172 GeV is used as input parameter. 


is given by A = when no cut is applied. The hnite- 
width infrared divergence of the non-resonant match¬ 
ing coefficient that cancels the corresponding ultravio¬ 
let divergence of the non-relativistic current correlation 
function appears as an endpoint divergence of the inte¬ 
gral above as pj inf approaches the on-shell value. 
At NLO, the divergence arises only from the diagram 
shown in Figure The integrand behaves as 


Hi 



Mw\ 

mf / 


const X 


1 


(mf-pf) 


,2i2’ 


(31) 


which leads to a divergent integral (301 in four dimen¬ 


sions. Dimensional regularization sets such linearly di¬ 
vergent integrals to hnite numbers, which explains the 
absence of explicit p„ scale dependence at this order. 

Similar to the case of VT-boson pair production, the 
leading non-resonant contribution to the top-pair cross 
section is a nearly energy-independent, negative correc¬ 
tion, which amounts to a few percent above the topo- 
nium peak, and to around 20% a few GeV below the 
peak. With an invariant-mass cut the NLO correction 
is shown as the lower (black) solid and dashed lines in 
Figure [T^ 


Since the largest sensitivity to the top-quark mass 
comes from the steep rise of the cross section below the 
peak, the non-resonant contributions at NLO and even 
NNLO are essential. The NNLO terms correspond to 
three-loop cut diagrams. The complete calculation has 
not yet been performed. However, in the presence of an 
invariant-mass cut satisfying F, «c AM, «c m,, the sin¬ 
gular terms as m,lAM, oo have been extracted in two 
different ways ll^ The calculation of llTSl . which 
starts from the non-resonant side, also conhrms explic¬ 
itly the cancellation of 1 je hnite-width divergence poles 
with the non-relativistic contributions. The upper set of 
lines (blue solid and dashed) in Figure [T^ shows that 
the NNLO correction is only about half as large than 
the NLO one. 

4. Summary and further results 

In this article I reviewed the treatment of unstable par¬ 
ticles in perturbative quantum held theory based on the 
scale hierarchy F M. Once scale separation is taken 
as the guiding principle, and an effective held theory 
is constructed, gauge-invariance and the consistency of 
the all-order resummation is automatic. The effective 
theory describes the scattering processes that leave the 
resonance close to its complex mass shell. Some aspects 
are therefore closely related to other effective theories 
that describe heavy particles close to their mass-shell. 
The concrete applications considered so far can be de¬ 
scribed by the master formula which captures res¬ 
onant production and decay, as well as all non-resonant 
“background” processes. 

The effective theory approach appears most fruitful, 
when it is applied to inclusive quantities, where it leads 
to particularly simple, even completely analytic results; 
to processes that require other resummations on top of 
the self-energy of the unstable particle; and to processes 
where high precision is required, beyond NLO accuracy, 
for which automated tools are not yet available. The 
line-shape and pair production near threshold discussed 
here are examples of such situations. 

As with other effective theories it is more difficult to 
predict differential distributions, unless the scales as¬ 
sociated with the observable can be assigned a unique 
scaling with respect to the small parameters that de- 
hne the EFT expansion. To circumvent this problem, 
a hybrid approach has been followed in If3^ [39l . ap¬ 
plicable to NLO calculations, in which the simplihca- 


^Further, the leading term in an expansion in the parameter p = 
1 - Mw/m, has been obtained in 03 and Gzl. with different results. 
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tions provided by the EFT are used in the virtual cor¬ 
rections (which have tree kinematics), while real emis¬ 
sion is computed in the full theory with the complex- 
mass prescription. A fully differential calculation has 
been done recently for the mixed (9(a'ema'i) corrections 
to Drell-Yan production ITOl employing a mixture of di¬ 
agrammatic and EFT-inspired techniques. 
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